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Abstract 

As  part  of  the  1  993  Wave  Basin  Experiment  (WABEX-93),  water  velocities  and 
gas  fluxes  were  measured  in  a  surf  pool  in  Irvine,  California,  under  plunging 
and  spilling  waves  with  heights  up  to  1 .2  m.  Vertical  profiles  of  turbulence- 
dissipation  rates  were  estimated  from  the  velocity  measurements  made  with 
an  Acoustic  Doppler  Velocimeter  (ADV).  The  ADV  measured  three  compo¬ 
nents  of  velocity  at  a  sampling  rate  of  25  Hz.  Turbulence-dissipation  rates  • 
were  estimated  from  the  magnitude  of  the  wavenumber  spectra  of  verti¬ 
cal  velocity  using  the  universal  form  for  the  inertial  subrange.  Dissipation 
rates  beneath  unbroken  waves  were  found  to  be  low  in  magnitude  and  uni¬ 
form  in  the  vertical.  Beneath  broken  waves,  dissipation  rates  were  found 
to  be  significantly  higher  in  magnitude  throughout  the  water  column  and 
to  increase  toward  the  water  surface.  Near-surface  dissipation  rates  be¬ 
neath  breaking  waves  were  found  to  follow  an  exponential  dependence  on 
whitecap  coverage  for  all  but  one  wave  height  condition.  This  suggests 
that  whitecap  coverage  can  be  used  to  scale  increases  in  dissipation  rate  for 
breaking  waves.  The  gas  transfer  velocity  due  to  turbulence  generated  by 
near-surface  currents  and  nonbreaking  waves,  and  the  gas  transfer  velocity 
due  to  turbulence  generated  by  breaking  waves  have  been  estimated  from 
their  respective  dissipation  rates.  These  transfer  velocities  have  been  used 
with  the  bubble-mediated  transfer  velocity  to  predict  the  total  gas  transfer 
velocity  in  the  surf  pool.  Predicted  gas  transfer  velocities  agree  well  with 
measured  gas  transfer  velocities. 


1  Introduction 

The  flux  of  gases  to  or  from  bodies  of  water  can  be  estimated  as  the  prod¬ 
uct  of  the  air -water  concentration  difference,  AC,  of  a  gas  and  its  transfer 
velocity,  ki.  Because  of  the  difficulties  associated  with  direct  oceanic  mea¬ 
surements  of  ki,  it  is  often  calculated  from  parameterizations  derived  em¬ 
pirically  from  wind  speed.  High  surface  winds  generate  surface  currents 
and  breaking  waves,  which  generate  bubbles  and  affect  mixing.  As  part  of 
an  experiment  designed  to  help  develop  a  method  for  estimating  ki  at  high 
wind  speed,  gas  fluxes  and  turbulence  were  measured  in  a  large  outdoor 
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siorf  pool.  A  more  complete  description  of  the  experiment  and  the  gas-flux 
measurements  can  be  found  in  Asher  et  al.  [1995]  and  Warminkhof  et  al. 
[1995]. 

Results  from  the  gas-exchange  measurements  suggest  that  ki  can  be 
partitioned  into  a  component  due  to  near-surface  turbulence  generated  by 
currents  and  nonbreaking  waves  (Um),  one  due  to  turbulence  generated  by 
breaking  waves  (kr)>  and  one  due  to  bubble-mediated  transfer  {ks).  For  gas 
evasion  with  AC  much  less  than  zero,  ki  can  be  written  in  terms  of  the 
fractional  area  of  whitecap  coverage,  Wc,  as 


ki  =  +  Wcikx  —  kw))  +  Wckh  (1) 

Asher  et  al.  [1995]  showed  that  (1)  could  be  applied  to  gas  transfer  mea¬ 
sured  in  the  surf  pool  using  an  exphcit  functional  form  for  kg  developed 
using  data  collected  in  a  whitecap  simulation  tank.  Although  (1)  provided 
reasonable  estimates  of  ki  in  the  surf  pool,  it  was  observed  to  systematically 
underpredict  ki  when  Wc  was  small.  Asher  et  al.  [1995]  hypothesized  that 
this  discrepancy  was  caused  by  incorrectly  parameterizing  the  dependence 
of  kw  and  kj  on  Wc- 

Equation  1  was  derived  by  assuming  that  k^  and  kj  were  constant  and 
not  functions  of  Wc-  Because  k^  and  kj  are  determined  by  the  vertical 
structure  and  mtensiry  of  the  near-surface  aqueous-phase  turbulence,  they 
are  constant  only  if  the  turbulence  in  the  surf  pool  is  constant  as  Wc  changes. 
However,  it  is  reasonable  to  expect  that  the  turbulence  associated  with  break¬ 
ing  waves  increases  with  wave  height  and  that  the  background  turbulence 
levels  in  the  surf  pool  wfll  increase  with  increasing  wave  energy.  The  present 
study  confirms  that  the  turbulence-dissipation  rate  increases  with  wave  height 
Furthermore,  the  increase  in  dissipation  rate  can  be  scaled  using  the  white- 
cap  coverage  of  the  breaking  waves,  and  inclusion  of  turbulence-dissipation 
rates  in  estimating  ku  and  kj  improves  the  ability  of  (1)  to  predict  gas  trans¬ 
fer  velocities. 


2  Experiment 

As  part  of  the  1993  Wave  Basin  Experiment  (WABEX-93),  estimates  of  turbu¬ 
lence-dissipation  rates  were  made  in  the  Hurricane  Harbor  surf  pool  at  Wild 
Rivers  Waterpark  in  Irvine,  California,  under  breaking  waves  with  heights  up 
to  1.2  m  (see  Asher  et  al.  [1995]  for  details  of  the  surf  pool  and  wave  gen¬ 
eration).  The  objective  of  the  measurements  described  here  was  to  obtain 
velocity  profiles  and  pressure  measurements  at  multiple  locations  in  the  surf 
pool  with  resolution  sufficient  to  estimate  profiles  of  turbulence-dissipation 
rate. 

Turbulence-dissipation  rates  were  estimated  from  velocity  measurements 
made  with  a  SonTek  acoustic-doppler  velocimeter  (ADV).  The  ADV  measured 
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Figure  1:  Plan  view  of  the  surf  pool  at  Wild  Rivers  Waterpark  in  Irvine,  California, 
showing  the  three  ADV  measurement  locations. 


three  components  of  velocity  in  a  0.5-cm^  sampling  volume  located  10  cm 
below  the  probe,  rninimizing  the  effect  of  the  sensor  on  the  flow  field.  Mea¬ 
surements  were  recorded  at  25  Hz  for  runs  lasting  approximately  300  s 
(70-100  waves).  Because  the  wave  field  was  highly  repeatable,  profQes  were 
made  by  compositing  data  from  three  to  seven  runs,  each  with  the  instru¬ 
ment  mounted  at  the  same  horizontal  location  but  at  a  different  elevation. 
Profiles  were  obtained  at  three  locations  across  the  surf  zone  under  spilling 
and  plimging  waves  generated  at  four  energy  levels  by  the  wavemaker,  with 
nominal  unbroken  wave  heights  of  0.3,  0.6,  0.9  and  1.2  m.  Figure  1  is  a 
schematic  drawing  of  Hurricane  Harbor  showing  the  three  ADV  measmre- 
ment  locations,  denoted  A,  B,  and  C.  Simultaneous  measurements  of  pres¬ 
sure  were  made  with  a  precision  strain-gage  located  on  the  bottom  directly 
beneath  the  ADV.  When  ADV  measurements  were  made  very  close  (<15  cm) 
to  the  bottom  the  pressure  transducer  was  located  ~  30  cm  to  the  side.  Video 
images  of  the  water  surface  were  recorded  and  later  used  to  determine  wave 
type  and  breaker  location. 


3  Data  Analysis 

Flow  beneath  waves  can  be  separated  into  three  components:  mean,  peri¬ 
odic  (wave),  and  fluctuatiag  (turbulence)  components.  The  objective  of  the 
data  analysis  was  to  quantify  the  turbulence  component.  Several  approaches 
were  considered.  In  laboratory  studies  with  highly  repeatable  waves,  aU  ve¬ 
locity  deviations  from  the  ensemble-mean  wave  velocity  can  be  attributed 
to  turbulence  (Nadaoka  et  al.  [1983]).  Although  waves  in  the  surf  pool  ap- 
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peared  qmte  regular,  the  detailed  pressure  measurements  showed  wave-to- 
wave  variations  in  pressure.  Irregularity  in  the  wave  field  was  most  likely 
caused  by  the  combined  effects  of  variability  in  the  wavemaker,  seiching  of 
the  pool,  reflected  waves,  and  episodic  outbreaks  of  a  rip  current.  Estimates 
of  turbulence  based  on  deviations  from  the  ensemble-mean  velocity  were 
precluded  because  they  would  overestimate  turbulence  by  including  non- 
turbulent,  pressure-induced  fluctuations.  An  alternative  approach  for  ex¬ 
tracting  turbulence  from  velocity  records  involves  the  portion  of  the  records 
associated  with  wave-induced  variations  of  the  water-surface  (Kitaigorodskii 
et  al.  [1983],  Agrawal  and  Aubrey  [1992]).  This  frequency-domain  method 
assumes  that  velocity  variations  that  are  coherent  with  pressure  fluctuations 
are  caused  by  two-dimensional  progressive  waves,  and  that  noncoherent  ve¬ 
locity  fluctuations  can  be  attributed  to  turbulence.  However,  the  technique 
is  prone  to  error  if  waves  have  directional  spread  (Herbers  and  Guza  [1993]). 
We  instead  used  the  inertial  dissipation  method  of  George  et  al.  [1994],  mod¬ 
ified  to  take  advantage  of  the  nearly  regular  wave  field  and  the  measurement 
capabilities  of  the  ADV. 


3.1  Data  Processing 

One  advantage  of  the  ADV  was  its  ability  to  rapidly  re-estabhsh  good  velocity 
measurements  after  the  loss  of  acoustic  signals  in  bubble  plumes  or  during 
subaerial  exposure  in  the  trough  of  a  wave.  Initial  processing  identified  and 
rejected  data  from  times  when  the  ADV  was  exposed  or  large  amounts  of 
bubbles  were  present.  A  prehminary  criterion  for  rejection  v/as  based  on 
the  correlation  coefficient  between  the  transmitted  acoustic  signal  and  the 
received  backscattered  signal.  Data-rejection  rates  increased  from  near  zero 
at  depth  to  as  much  as  37  %  near  the  siuface  under  larger  waves.  Acceptable 
data  for  70  to  100  successive  waves  were  phase  aligned  using  the  peak  in  the 
cross  covariance  of  the  slope  of  the  pressure  signal,  and  ensemble  means 
and  standard  deviations  for  three  velocity  components  and  pressure  were 
estimated  (Figure  2).  An  iterative  check  on  data  was  then  performed:  data 
were  rejected  if  the  velocity  magnitude  exceeded  3.7  standard  deviations 
about  the  ensemble  mean  at  each  phase  in  the  wave,  and  ensemble  statistics 
were  recalculated.  Fmahy,  mean  and  root-mean-square  (RMS)  velocities  were 
determined  for  acceptable  data  from  each  data  run. 


3.2  Estimates  of  Dissipation  Rates 

In  steady  flow  with  isotropic,  fcdly-developed  turbulence,  kinetic  energy  is 
transferred  from  the  mean  flow  to  large  eddies,  then  to  small  eddies,  and 
is  finally  dissipated  by  viscosity.  Under  these  conditions,  the  turbulence- 
dissipation  rate  can  be  estimated  by  the  magnitude  of  the  wavenumber  spec¬ 
tra  in  the  inertial  subrange,  which  takes  the  form: 
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E(K)  =  (2) 

where  k  is  wavenumber,  E(k)  is  wavenumber  spectral  density,  a  is  the 
turbulence  dissipation  rate,  and  a  is  the  Kolmogorov  constant.  The  inertial 
subrange  extends  from  large  eddies,  the  scale  of  which  is  typically  deter¬ 
mined  by  physical  dimensions  of  the  flow  (e.  g.,  depth)  to  the  Kolmogorov 
microscale,  which  is  determined  by  kinematic  viscosity  and  dissipation  rate. 

Estimates  of  dissipation  rate  were  made  for  the  surf  pool  data  using  (2), 
and  assume  that  the  turbulence  was  fully  developed  and  isotropic.  The  ver¬ 
tical  component  of  the  ADV  velocity  signal  was  used  to  calculate  dissipation 
rates  because  it  had  the  highest  signal-to-noise  ratio  of  the  three  compo¬ 
nents.  Spectral  estimates  were  obtained  using  a  windowing  and  ensemble- 
averaging  (Welch)  method.  Because  data  were  removed  during  initial  pro¬ 
cessing,  the  spectra  for  each  windowed  time  series  was  computed  using  a 
method  appropriate  for  irregularly  spaced  data  (Press  et  al.  [1992]).  Entire 
data  windows  were  omitted  if  more  than  50%  of  the  data  points  were  miss¬ 
ing.  Spectral  estimates  were  obtained  by  first  detrending  the  data  by  sub¬ 
tracting  out  the  ensemble-mean  periodic  (wave)  velocity  component,  then 
computing  spectra  on  a  windowed  data  segment  extending  over  three  suc¬ 
cessive  waves,  and  finally  ensefnble-averaging  the  spectra.  Windowed  seg¬ 
ments  were  overlapped  67  %.  This  procedure  yielded  spectral  estimates  with 
resolved  bandwidths  of  0.16  Hz  and  156  degrees  of  freedom  for  typical  time 
series  containing  80  waves. 

The  time-series  data  at  a  fixed  point  were  used  to  estimate  the  energy- 
-  density  spectrum  as  a  function  of  frequency,  not  wavenumber.  The  fre¬ 
quency  spectra  were  converted  to  wavenumber  spectra  using  Taylor’s  hy¬ 
pothesis  of  frozen  turbulence,  assuming  that  the  time  scale  of  turbulent 
fluctuations  was  long  compared  with  the  time  scale  of  the  motion  advecting 
the  eddy  past  the  sampling  point.  The  RMS  orbital  velocity  was  used  as  the 
advective  velocity  for  the  conversion  from  measured  frequency  spectra  to 
wavenumber  spectra  (.Agrawal  et  al.  [1992]).  An  example  of  the  wavenumber 
spectra  generated  from  the  vertical  velocity  data  is  shown  in  Figure  3. 

Turbulence-dissipation  rates  (a)  were  estimated  from  the  magnitude  of 
the  wavenumber  spectra  in  the  inertial  subrange  using  (2).  For  each  data  run, 
the  wavenumber  spectra  were  individually  inspected  and  the  magnitude  of 
the  spectra  in  the  range  over  which  the  theoretical  -5/3  slope  was  fomd  was 
used  to  estimate  the  dissipation  rate  (typically  between  radian  wavenumbers 
of  0.4  and  1.5  cm"^). 

4  Turbulence  Dissipation  Results 

Profiles  of  a  beneath  plimging  waves  are  shown  in  Figures  4a  and  4b,  and 
beneath  spilling  waves  in  4c.  Elevation  is  normalized  by  the  water  depth. 
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Time  (s) 


Figure  2:  Ensemble-average  velocities  and  pressure  under  plunging  wave  of  nominal 
height  0.9  m  at  location  A,  53  cmab:  a  Alongshore  velocity,  b  Cross-shore  velocity,  c 
Vertical  velocity,  d  Pressure. 


so  that  a  normalized  elevation  of  1.0  is  at  the  mean  water  surface.  Under 
unbroken  waves,  a  profiles  were  generally  low  in  magnitude  and  vertically 
uniform.  Beneath  plunging  waves  measured  at  location  A  (Figure  4a),  a  was 
highest  near  the  surface  and  increased  with  increasing  wave  height.  Dissipa¬ 
tion  rates  were  also  seen  to  increase  vdth  increasing  wave  height  throughout 
the  water  column. 

Profiles  of  a  at  location  B  (Figure  4b)  are  similar  to  profiles  at  location  A 
(Figure  4a),  although  the  mid-water  column  a  are  higher  at  B,  and  there  is  at 
B  little  difference  in  a  between  the  0.9-m  and  1.2-m  wave  conditions.  Strong 
rip  currents  at  B  could  be  responsible  for  the  increase  in  the  a  for  the  mid¬ 
water  column  compared  with  A.  Alternatively,  the  increased  distance  from 
breaking  (cross-shore  location)  at  B  could  have  allowed  homogenization  of 
turbulence  through  the  water  column.  Similar  values  of  a  for  wave  heights 
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Wavenumber,  1/cm 

Figure  3:  Representative  wavenumber  spectrum  of  vertical  velocity  showing  inertial 
subrange  slope)  under  plunging  wave  of  nominal  height  0.9  m  at  location  A, 

53  cmab. 


of  0.9  m  and  1.2  m  could  indicate  a  saturation  value  for  f  at  some  distance 
from  the  break  point.  At  B,  velocity  measurements  were  made  very  close 
to  the  pool  bottom,  and  observed  increases  in  f  are  likely  associated  vvdth 
shear-induced  turbulence  in  the  bottom  boundary  layer  (Figure  4b). 

Beneath  spilling  waves  at  location  C  (Figure  4c),  profiles  of  f  show  more 
variation,  possibly  because  spilling  waves  were  less  regular.  Like  e  measured 
beneath  plunging  waves  at  B,  e  was  high  in  the  middle  of  the  water  column 
but,  in  contrast  to  measurements  of  £  at  B,  f  measured  beneath  spilling 
waves  at  C  did  not  appear  to  reach  satuiction  near  the  surface. 


4.1  Characteristic  Dissipation  Rates 

Dissipation  rates  from  two  depths  were  used  to  characterize  turbulence  gen¬ 
erated  by  breaking  waves  (et),  and  background  mechanical  turbulence  gen¬ 
erated  by  currents  and  nonbreaking  waves  (em)-  The  Et  was  taken  to  be  the 
E  interpolated  at  a  normalized  elevation  of  0.8  for  the  plunging  waves,  and 
0.65  for  the  spilling  waves  (upper  limit  of  measurements).  These  depths  are 
consistent  with  the  assumption  that  the  dominant  length  scale  for  turbu¬ 
lence  generated  by  wave  breaking  is  proportional  to  the  breaking  wave  am¬ 
plitude.  It  was  assumed  that  the  characteristic  ej  and  Em  for  gas  exchange  in 
the  entire  surf  pool  were  less  than  et  and  em  at  location  A  and  greater  than 
Et  and  em  at  location  B  beneath  plunging  waves,  and  that  the  characteristic 
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Figure  4:  Vertical  profiles  of  turbulence-dissipation  rate,  e,  beneath  a  plunging  waves 
at  location  A,  b  plunging  waves  at  location  B,  and  c  spilling  waves  at  location  C,  of 
nominal  waveheights  0.3,  0.6,  0.9,  1.2  m. 


et  and  em  were  approximately  equal  to  ej  and  Em  at  location  C  beneath  the 
spilling  waves.  Values  for  ej  are  listed  in  Table  1  for  che  two  wave  types  and 
four  wave  heights.  The  dissipation  rate  for  turbulence  generated  by  non¬ 
breaking  waves  and  currents,  Em^  was  determined  at  a  normalized  elevation 
of  0.2  (Table  1).  This  level  was  generally  below  the  region  of  large  gradi¬ 
ents  in  the  e  profiles.  Turbulence  generated  in  the  bottom  boundary  layer 
was  observed  below  this  level  (Figure  4b),  but  does  not  appear  to  influence 
E  above  normalized  elevations  of  0.1. 


4.2  Dissipation  Rates  and  Whitecap  Coverage 

■  The  area-weighted  fractional  coverage  of  actively  breaking  waves,  Wq  ,  was 
measured  in  the  pool  (Asher  et  al.  [1995]),  and  values  are  listed  in  Table 
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Table  1:  Turbulence-dissipation  rates  and  whitecap  coverage 


Wave  T  /pe 

Nominal 

Height 

£m  (cmVs^) 

£t  (cmVs^) 

VKc  (%) 

Plunging 

1.2  m 

75 

499 

7.15 

Plunging 

0.9  m 

45 

299 

3.96 

Plunging 

0.6  m 

2 

66 

2.01 

Plunging 

0.3  m 

0.6 

0.8 

0.00 

Spilling 

1 .2  m 

65 

428 

2.90 

Spilling 

0.9  m 

21 

152 

2.51 

Spilling 

0.6  m 

5 

6 

0.54 

Spilling 

0.3  m 

1 

2 

0.00 

1.  Figure  5  is  a  log-normal  plot  of  fr  versus  Wc  from  the  data  in  Table 
1.  £  for  spilling  waves  at  all  four  heights  and  ej  for  breaMng  waves  with 
heights  of  0.3  m,  0.6  m,  and  0.9  m  were  found  to  follow  the  same  exponential 
dependence  on  Wc-  This  suggests  that  Wq  can  be  used  to  scale  increases  in 
£j  caused  by  breaJdng  waves.  It  is  not  clear  why  fr  for  the  1.2-m  plunging 
breaking  waves  does  not  follow  the  same  dependence  as  the  other  seven 
wave  conditions,  but  it  should  be  noted  that  these  conditions  presented 
the  most  difficult  data-analysis  problems,  and  many  data  points  during  the 
times  of  high  velocities  were  eliminated.  As  a  result,  maximum  dissipation 
rates  could  have  been  underestimated. 

5  Parameterization  of  Gas  Transfer  Velocities 

Lamont  and  Scott  [1970]  showed  that  for  gas  exchange  across  an  unbroken 
water  surface,  ki  could  be  estimated  using  £.  Using  their  relation,  ku  and 
kj  can  be  written  as 


kx  =BSc-^i^{£xv)^'^  {X  =  M,T)  (3) 

.  where  B  is  a  dimensionless  constant.  Sc  is  the  Schmidt  number  of  the  gas 
(equal  to  the  ratio  of  kinematic  viscosity  to  molecular  diffusivity),  v  is  the 
kinematic  viscosity  of  water,  and  sm  and  fr  are  the  dissipation  rates  of  the 
wave-current-generated  turbulence  and  the  breaking-generated  turbulence, 
respectively.  In  support  of  (3),  Asher  and  Pankow  [1986]  foimd  that  ki  for 
transfer  of  carbon  dioxide  (CO2)  through  a  cleaned  water  surface  scaled  lin¬ 
early  with  They  also  found  that  B  was  equal  to  0.65. 

Asher  et  al.  [1995]  measured  ki  for  CO2,  helium  (He),  nitrous  oxide  (N2O), 
and  sulfur  hexafluoride  (SFe)  as  a  function  of  Wc  in  the  surf  pool.  They 
developed  an  empirical  parameterization  of  ki  using  (1)  by  assuming  kM  and 
kr  were  constant  wath  increasing  Wc  and  taking  ks  =  As 
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Figure  5:  Whitecap  coverage,  Wq,  versus  surface  turbulence-dissipation  rate,  sr,  for 
plunging  waves  at  locations  A  and  B  and  spilling  waves  at  location  C. 


can  be  seen  in  Figure  5  however,  sm  and  et  were  not  constant  wdth  increasing 
Wc,  and  (3)  shows  fcw  and  kj  vsill  also  increase.  Therefore,  it  is  not  surprising 
that  Asher  et  al.  [1995]  found  that  parameterizing  ki  using  (1)  and  assuming 
that  kjyf  and  kr  were  constant  was  unsatisfactory. 

The  dissipation  rates  listed  in  Table  1  were  used  to  calculate  kn  and  fcj 
using  (3)  and  assuming  B  would  be  a  constant  fit  by  the  data.  These  transfer 
velocities  were  substituted  into  (1)  to  account  for  the  effect  of  the  increase 
in  turbulence  with  increasing  Wc-  Written  explicitly,  ki  in  the  wave  basin 
can  be  parameterized  as 


ki  =  +  Wc  ((£rv)^^^  -  +  biWcOi  (4) 


The  resulting  model  was  then  fit  using  least-squares  linear  regression  to 
the  measured  ki  data  from  Asher  et  al.  [1995]  to  determine  the  coefficients 
bi  and  B  in  (4).  This  fit  results  in  B  =  0.58  and  bi  =  6.30  m  s”^  with  a 
coefficient  of  determination  of  0.84.  The  value  of  B  found  for  the  WABEX- 
93  data  is  in  good  agreement  with  that  determined  by  Asher  and  Pankow 
[1986]. 

Figure  6  shows  the  result  of  fitting  (4)  to  the  surf  pool  data  by  plotting 
ki  for  CO2,  He,  N2O,  and  SFe  calculated  using  (4)  versus  the  experimentally 
determined  value.  There  is  excellent  agreement  between  the  calculated  and 
measured  ki  values.  Asher  et  al,  [1995]  found  a  reduced  chi-squared  (x^)  of 
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192  for  the  constant  and  kj  parameterization.  The  for  the  data  in  Fig¬ 
ure  6  was  70,  which  shows  that  including  the  effects  of  increasing  turbulence 
with  increasing  Wc  improves  the  predictive  accuracy  of  a  parameterization 
based  on  (1). 

The  parameterization  given  in  (4)  can  also  be  used  to  estimate  the  fraction 
of  ki  that  is  due  to  turbulence  and  the  fraction  of  ki  that  results  from  bubble- 
mediated  processes.  For  SFe  for  a  0.9-m  spiOing  wave,  turbulence  processes 
accounted  for  38%  of  the  total  transfer  velocity  and  bubble  processes  were 
62%  of  the  total.  In  contrast,  the  turbulence  fraction  for  CO2  imder  the 
same  conditions  was  46%  and  the  bubble  fraction  was  54%.  The  decreased 
importance  of  bubble-mediated  transfer  for  CO2  compared  with  SFe  agrees 
with  the  model  predictions  of  Memery  and  Merlivat  [1985].  Their  results 
indicated  that  ks  for  a  soluble  gas  such  as  CO2  would  be  less  than  ke  for  an 
insoluble  gas  like  SFe. 

6  Conclusion 

Turbulence  dissipation  rates  were  found  to  increase  with  increasing  wave 
height  for  spilling  and  plimging  waves.  The  surface  dissipation  rates  scaled 
as  an  exponential  fimction  of  whitecap  coverage  for  spilling  waves  with  nom¬ 
inal  heights  of  0.3  m,  0.6  m,  0.9  m,  and  1.2  m  and  for  plunging  waves  with 
nominal  heights  of  0.3  m.  0.6  m  and  0.9  m.  For  reasons  that  are  not  clear  at 
present,  the  surface  dissipation  rate  for  plunging  waves  with  a  height  of  1.2 
m  did  not  follow'  this  exponential  scaling. 

Asher  et  al.  [1995]  found  that  parameterizing  ki  for  CO2,  He,  N2O,  and  SFe 
assuming  ky,  and  kj  were  constants  did  not  provide  a  satisfactory  fit  with 
the  measured  surf  pool  data  =  192).  Substitution  of  the  direct  estimates 
of  ky  and  kj  from  (3)  and  dissipation  rates  estimated  from  the  surf  pool 
measurements  improved  the  agreement  between  measured  and  predicted 
ki  values  =  70). 

In  the  surf  pool,  increases  in  Wc  were  caused  by  increases  in  height  of 
the  breaking  wave.  It  could  be  expected  that  a  similar  process  occxirs  under 
oceanic  conditions,  where  increasing  wind  speed  leads  to  larger  wave  size 
and  higher  whitecap  coverage.  The  success  of  the  approach  described  here 
suggests  that  successful  application  of  (1)  to  ocean  conditions  will  require 
knowledge  of  ks  and  the  dependence  of  ky  and  kj  on  wind  speed.  Although 
this  comphcates  the  parameterization  of  ki  in  the  ocean,  it  will  provide  an 
accurate  estimation  of  the  effect  of  the  turbulence  and  bubbles  generated 
by  breaking  waves  on  air-sea  gas  fluxes. 
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